Optical conductivity of unconventional charge density wave systems: 

Role of vertex corrections 
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The optical conductivity of a d-CDW conductor is calculated for electrons on a square lattice and 
a nearest-neighbor charge-charge interaction using the lowest-order conserving approximation. The 
spectral properties of the Drude-like peak at low frequencies and the broad hump due to transitions 
across the gap at large frequencies are discussed, also as a function of temperature and of the second- 
nearest neighbor hopping term t' . We find that vertex corrections enhance the d.c. conductivity, 
make the Drude peak narrower and provide a smooth transition from a renormalized regime at 
low to the bare theory at high frequencies. It is also shown that vertex corrections enhance the 
temperature dependence of the restricted optical sum leading to a non-negligible violation of the 
sum rule in the d-CDW state. 

PACS numbers: 71.45.Lr,78.20.Bh,74.72.-h 



I. INTRODUCTION 

The interaction between electrons may give rise to 
three different types of order parameters describing 
superconductivity, magnetism or a modulation of the 
charge density. As an example let us consider the Heisen- 
berg interaction, 



J 



t 
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The sum in Eq. (JQ runs over nearest neighbor sites i and 
j, u is the vector of Pauli matrices, and c\ a , Ci a are cre- 
ation and annihilation operators for electrons at site i 
with spin projection a. In a mean-held treatment the 
expectation value of the two creation (or annihilation) 
operators, i.e., the order parameter, may be non-zero 
which leads to singlet superconductivity with extended 
s-wave or d-wave symmetry or triplet superconductivity 
with p-wave symmetry. If the average of two operators 
from the same site is non-zero we have a modulation of 
finite local spin moments describing a magnetic state. 
If the expectation value of two operators with the same 
spin projection is non-zero we deal with a charge-density 
wave (CDW) with an internal symmetry, because the 
two operators in the order parameter belong to different 
sites. This case will be considered in the following. For 
a square lattice with 4 nearest neighbors there are four 
different local order parameters associated with each of 
the four bonds which are in general complex valued but 
also subject to the condition that the resulting mean-field 
Hamiltonian must be Hermitian. The imaginary parts 
correspond to circulating currents expressing the break- 
ing of time reversal symmetry»i The order parameters of 
a CDW may be classified by the point group symmetry 
and have in our case extended s-wave, d-wave or p-wave 
symmetries. Moreover, they show a modulation through- 
out the crystal with a wave vector Q corresponding to 



the total momentum of electron-hole pairs. CDW sys- 
tems with an internal structure are called unconventional 
charge-density waves^ to distinguish them from conven- 
tional ones where the electron and hole sit at the same 
site and do not exhibit any internal structure. 

Unconventional charge-density wave systems have 
been invoked to explain experimental features in organic 
conductors^ 2D transition metal dichalcogenides^ and 
high-T c superconductorsi£i£iL2i2ii& In particular, it has 
been shown that a charge density wave with d-wave sym- 
metry (d-CDW or, shortly, DDW) competes with d-wave 
superconductivity in the t — J model in the large- A limit 
where N is the number of spin projections. 6 The phase di- 
agram of this model has a quantum critical point (QCP) 
separating at T = the normal phase at large doping 
from a d-CDW state at lower doping if superconductiv- 
ity is disregarded. Allowing also for superconductivity 
the QCP separates a pure superconducting state from 
a ground state where superconductivity and a d-CDW 
coexist. Due to the presence of the d-CDW the transi- 
tion temperature for superconductivity exhibits a max- 
imum, to be identified with optimal doping, and then 
decreases strongly in the underdoped regime. The d- 
CDW state exists above T c near optimal and the un- 
derdoped region which means that the pseudogap phase 
should be identified with the d-CDW in this model. Var- 
ious observable quantities have been calculated within 
this picturo 1 ! 6 ! 8 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 and compared with 
experimental data but no final conclusions have been 
reached. 

In the following we will study the optical conductivity 
of a d-CDW system. In contrast to other recent workiiiiS 
vertex corrections will be taken into account within a con- 
serving approximation schemed Otherwise a proper def- 
inition of the current already poses a problem. Within a 
mean-field theory the current can be defined via the con- 
tinuity equation for the charge density. The explicit cal- 
culation of the time derivative of the charge density using 
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the mean-field Hamiltonian of a d-CDW yields the follow- 
ing expression for the longitudinal current in k spacejSi 

j( q ' f ) = lu H [( V kek)4_ q/2cr c k+q/2CT 

kcr 



-i(VkA k )4_ q/2CT c k+ Q +q /2 



(2) 



e k and A k are one-particle energies and the d-wave or- 
der parameter, respectively. The first term in Eq.Q is 
the standard expression for the current, the second one 
is unconventional and due to the momentum dependence 
of the order parameter. The use of the current definition 
Eq.|2j| in calculating the optical conductivity guarantees 
that the continuity equation is satisfied 21 but it is not 
apriori clear whether it is at finite frequencies more cor- 
rect than the usual procedure based only on the first 
term in Eq.(J2Jl. In particular, one can argue that the ex- 
act interaction term in the Hamiltonian should be used 
in defining the current. One would then obtain only the 
first term in Eq.(j2J because the interaction, for instance, 
a Heisenberg term, commutes with the density operator. 
Below we will show that a correct treatment defines the 
current only via the first term in Eq. J2J) but takes vertex 
corrections (VC) into account. In this way the effective 
current at low-frequency is of the form of Eq. © and sat- 
isfies the continuity equation whereas it resembles more 
a bare current at high frequencies. 

The paper is organized as follows. In the second sec- 
tion we introduce some properties of the Hamiltonian of 
a d-CDW system. In section III we derive and solve the 
equation for the current vertex. In section IV the opti- 
cal conductivity is calculated and several limiting cases 
such as high and low frequencies are considered as well 
as the weight of the Drude peak, optical transitions and 
the optical sum rule. Section V contains results from a 
numerical evaluation of the obtained expressions, and we 
present our conclusions there. 



II. HAMILTONIAN 

We consider the following Hamiltonian for electrons on 
a square lattice, 



krr 



\ J ( k - k ')cL c k'<T'C k , v ,c k _ k / +k " (T . (3) 



kk'k" 
aa 

°La an< ^ Ck<T are creation and annihilation operators for 
electrons with momentum k and spin projection <r, re- 
spectively. £ k °^ is equal to e k — /i, where e k are the bare 
one-particle energies and /i the chemical potential. e k is 
given by e k = — 2t(cos(k x ) +cos(fc y )) + 4t'cos(k x )cos(k y ), 
where t and t' denote hopping amplitudes between near- 
est and next-nearest neighbors, respectively. The second 



term in TL is a Heisenberg interaction between nearest 
neighbors with J(k) = 2J(cos(k x ) + cos(fc a )), measuring 
lengths in units of the lattice constant. From now on we 
also will measure energy in units of 2 J which will simplify 
subsequent equations considerably. 

Extending the spin projections from two to N the 
ground state of the above Hamiltonian is a charge-density 
wave with d-wave symmetry for some doping regime away 
from half-filling and large TV. The corresponding wave 
vector is approximately Q — {tt,tt). The DDW order 
parameter is given by 

A k = A 7d (k) = -2 l7d (k)]T J dk' 7d (k')<4, + Q ff <c k v<}, 

a' 

(4) 

with 7ti(k) = (cos^) — cos(/c y ))/2. Here and below our 
measure of integration is 



dk — dk x dk y / {2ir) 2 



(5) 



In the following we assume that the parameters are 
such that this DDW state is realized and that only 
the charge channel is relevant in accordance with the 
large N limit. It is convenient to introduce the spinors 
\& kCT = (c kcr , c k+ Q CT ). The mean-field part of the Hamil- 
tonian can then be written as 



Hmf = ^'^L-Mka, 



with 



k.tr 



6c «A k 
-iA k £ k+Q 



(G) 



(7) 



The prime at the summation sign means that the sum 
extends only over the magnetic zone which is 1/2 of the 
original Brillouin zone. The mean-field contribution of 
the interaction term in Eq.Q contains a diagonal part, 

which renormalizes £ k °^ and yields £ k , and a non-diagonal 
part, described by A k . For the following it is convenient 
to introduce the abbreviations 



and 



e± = (& ± £k +Q )A 



£ k =[^+A k ] 1/2 . (9) 
-ff k can be diagonalizcd by the unitary transformation U, 

U = cxp(icr 1 6' k /2), (10) 
where # k is uniquely defined by the equations 

cos0 k = £-/£k, sin# k = A k /£ k . (11) 

a 1 denote Pauli matrices. We have UHU< = 
diag{e\, £2) = h and the new quasiparticle energies are 



(8) 



£1,2 — (,+ ± Ek- 



(12) 
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In the following we will deal will non-magnetic states 
which means that we may drop the spin label a and ac- 
count for sums over spin by simply inserting a factor 2. 
The fermionic Green's function matrix is given by 

G k (iu; n ) = (iu; n -H)-\ (13) 

where uj n denotes the (fermionic) Matsubara frequencies 
(2n + 1)ttT. Gk is diagonalized by the same matrix U, 
so that G — U^gU with g = (iui n — 

III. VERTEX EQUATION 

The integral equation for the current vertex is depicted 
graphically in Fig. ^where the two solid lines with arrows 
stand for the Green's function matrix G-\Jj,LO n ). Analyti- 
cally, the integral equation reads, 

f Q (q,zc n ) = f(q)+T^ f dkj(k-q) 

xi J 

x G k (M ; )f Q (k,iw„)G k (ia; ; + icu n ), (14) 

where the current vertex r Q (q, iu) n ) is a 2x2 matrix with 
the elements rg-(q, iu) n ). The bare current vertex 7 a (q) 



is diagonal and given by 7 a (q) = diag{d a ^i L , #a£ q +Q)- 
The integration over k runs over the chemical Brillouin 




FIG. 1: Diagrammatic representation of the integral equation 
for the current vertex 

zone, while q lies within the magnetic Brillouin zone. By 
choosing this convention we take into account all allowed 
scattering processes. Ea. lfHE}) can also be obtained by 
rewriting the interaction term in Nambu notation, using 
the diagrammatic rules and then extending the integra- 
tion over k throughout the chemical and not only the 
magnetic Brillouin zone. In Ea. l|14l) iuj n = 2ninT is the 
(bosonic) external frequency and the external momentum 
has been put to zero. We neglected retardation effects on 
the interaction amplitude J(k) which implies that I™ de- 
pends only on one frequency variable. 

Performing an analytical continuation Ea. (|14|l be- 
comes, 



r a (q, W ) = 7 Q (q)- 



dx r / •*"'- 

— dk J(k-q) n F (x) (G^{x)T a (k,oj)G^(x + uj) 

2m L V 



G%{x - w)r Q (k, u)Gg(xj) + (n F (x + u) - n F (x))G£(x)r a (k, lu)G£(x + lu) 



(15) 



Gj^ and G^ are the advanced and retarded Green's function matrices, respectively, and np the Fermi function. First 
we consider the static case, u = 0, where the last term vanishes. Inserting V£iJk for the vertex we obtain for the 
integral on the right-hand side of Ea. (fTHjl . 



dx 
2tti 



-JkJ(k - q)n F (x) (G k 1 (x)V^ k G k 1 (x) - G^V^kG^x)) 



(16) 



dx 



— dkJ(k-q)n,(a;)V£ (G£(x) - G£(x) } = V" / dxdk J(k - q)n F (x)G£(x) = -V°E q , (17) 



where 27riG k (x) = G£(x) - Gg(x). Noting that the 
mean- field Hamiltonian is V q ff q = 7 a (q) + V q £ q , we 
observe that (|15|1 has the solution^! 

f Q ( q ,o) = v°i/ q . (is) 

EcL ljlSJI is a Ward identity relating the vertex and the 
self-energy in the low-frequency limit. It is a consequence 
of the continuity equation for charge and current densi- 
ties. Since the self-energy matrix is always non-zero in 
the DDW state it means that the current vertex is always 
renormalized and never equal to the bare one. Only then 



the continuity equation can be fulfilled at low frequen- 
cies. This directly shows the importance of VC in the 
DDW state. However, there may be special cases of in- 
teractions where VC vanish. Such a case has recently 
been discussed^ where J was assumed to have the sep- 
arable form J(k, q) ~ /(q)/(k), with an even function 
/(q) = /(— q). For this choice of interaction the inte- 
grand in (|15JI is odd in k, hence the integral vanishes and 
VC are absent. 

We proceed now to solve the vertex equation Ea. i|14|) 
in the general case. The second term on the right-hand 
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side of l|14(l becomes after an analytic continuation, 



dkJ(q — k) / dx\dxi 



n F {x\) - n F (x 2 ) 
uj — X\ + X2 + iO 



xG^(xi)r a (k,w)G^(x 2 ). 



(19) 



the Green's function. The Green's functions take then 
the form g R (x) = diag[g^ k (x), gjf^ix)], where gf k {x) = 
(x — £j,k + */(2t)) _1 . A standard calculation gives for 
the 11 element of / in the limit of small r _1 <C T, 



Using the unitary transformation U of Eg. (|1U|I the 
Green's function matrix can be written as Gk = U^g^U^ 
with a diagonal t/k- The product of the three matrices 
in (|T§|) can be represented as U^g"{xi)T a g"{x2)U with 
fa _ uf^uK Next we expand f Q in terms of Pauli 
matrices 



pa 



(20) 



, r ,, r ,,, / , , n F (xi) - n F {x 2 ) „ , , „ . . 
( r + r s) / dxidx 2 - ; - ; Ar . gi^{x l )g 1M {x 2 ) 



' UJ — X\ + x 2 + iO ' 



n F (£l,k) > r a xSa \ 
"I 1 t 1 3 /' 



1 — iCJT 



(23) 



where n^, (x) is the derivative of the Fermi function. The 
non-diagonal components of I are given by the expression 



The s— wave component of the current, T" , is in general 
allowed by Eq. (|15fl but vanishes if J is a nearest- neighbor 
interaction. Hence we restrict ourselves to three compo- 
nents and write T a = Tq en + T^a^ + 02, and similarly 
for the transformed vertex, T a = r^oo + ^3^3 + f 1 "^. 
Using the explicit form for U we obtain 



pa 
1 2 



cos ffc — sm t>k 
sin 6>fc cos 0fe 



T 



(21) 



The diagonal components of the expression 



1= I dx ld X2 nF{Xl) : nF{X2 l r(^g''(x2) (22) 



UJ 



x\ + x 2 + iO' 



are regularized at uj — > in the presence of impurity 
scattering. We adopt the simplest picture of point-like 
scatterers, characterized by the inverse quasiparticlc life- 
time r. Note that isotropic scattering from such impu- 
rities allows us to neglect VC and to include t only in 



a 2- 



pa 
1 2 



»f(£i,k) ~ n_F(e 2 ,k) 
uj + i/r - £i,k + £2,k 



(1^2) 



Introducing shorthand notations, 



z 

Pi 
P± 

K 

we have 



{1-iujt)-\ 

-n' F (e jM ), 

(Pi±P 2 )/2, 

n F (ei^) ~ n F (e 2 ,u) 
(£2,k - £i,k) 2 - (w + i/ r ) ; 



(24) 
(25) 
(26) 

:(£2,k-£i,k), (27) 



/ = (To2(p + r« + P _r«) + ds^P-r" + P+ r«) + <7 2 #r£. 

(28) 

The next step in the evaluation of 119|) is to carry 
out the rotation U'lU. Suppressing the dependence on 
the external frequency in the notation we obtain from 
Eci. (|14|) . the following system of coupled equations, 



io(q) = 7o(q) + * / dkj(q-k)( P+ r^ + p_r«) 



r?(q) 
r?(q) 



7 3 "(q) 




y dkJ(q-k)f f 



r 



«P-r?(k) 




z P+ owrf(k) 
x J lr§(k) 



(29) 
(30) 



A further simplification to Eqs. I|29|) . (|30|l results from 
the nearest-neighbor form of the interaction J(k). In this 
case the integral in l|29[) vanishes, which can be seen by 
considering a shift q — > q + Q ; upon this shift 7q is un- 
changed and J(k — q) — * — J(k — q). The "symmetrical" 
component of the current is thus not renormalizcd and 
givenbyr£(q)=7«(q) = V«£+- 

In order to solve Eq. I|30|) we note that 73 (q) transforms 
like a vector, i.e., according to the representation T5 of 
the point group C^ v . Decomposing J(k — q) into sepa- 



rable kernels only sin(fc a ) sin(g Q ) can contribute because 
sin(fc a ) is the only function appearing in the decompo- 
sition which transforms according to the representation 
T$ . Taking into account the transformation properties of 
T under Ci v one finds that this argument applies also to 
T%. The solution of Ea. (|3llj) has thus the form, 



rf(q) = 73(q)+sm(<ZaK, 
r£(q) = sin(ff a )t§. 



(31) 
(32) 



The coefficients 



-3 ' u 2 



depend still on the external fre- 
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quency but no longer on momentum. It is convenient to 
introduce the notation 



7o 



P-7o 




73 



7V 



zp + 



v = tU Y ^ 1 r. 



and 



dk sin(fc Q )(zf t 7o + T>73), 
dk sin 2 (fc Q )U. 



Using this notation and the vectors 

d 



3 

a I ; 

2 , 



one represents Eq. (|30H in the following form 

c = d + Xc, 



(33) 
(34) 

(35) 
(36) 

(37) 



Using Ea. (|llf> more explicit expressions for the ele- 
ments of V and for d can easily be obtained, 



V U = (zp+£ + KAi)/El 

V 12 = V 21 =^A k (K-z P+ )/El 



U 



22 



(zp + Al + I<e)/El 



and 



4 = 



dksm(/e Q ) 
dk sin (k a ) 



^ h 7 3 Kn 



, P-7o"A k 



+ 73 ^2! 



IV. OPTICAL CONDUCTIVITY 



(39a) 
(39b) 
(39c) 



(40a) 
(40b) 



The frequency-dependent conductivity is given by one 
current loop with one renormalized and one bare vertex. 
(38) Analytically it is given by the following expression, 

I 



a{u)=uj- l Im [dk [ dxidx 2 ^^ ^^4rr[ 7 a (k)G^(a;i)f a (k,w)G^(a; 2 )] 

J J uj - xi + x 2 + «0 

Eq. H41|) allows the same analysis as was done above for the vertex, (cf. (|19|l ). The integration over xi, 2 yields 



(7{L0 = UJ 



- l Im J dkTr[j a (k)U^IU}, 



(41) 



(42) 



with / given by 128(1 . The integration is over the magnetic Brillouin zone. The trace in the last formula is represented 
as 



z^(k)(p +1 S + P-T%) + (7 3 a (k),0)f t 



zp_ 7o «(k) 




z P+ oWr§<(k) 
o k lrf(k) 



(43) 



Eas. (|42|) and (|43|l can be recast in the form 

a(uo) = ImQ/uo, (44) 
Q = z(a + 2a 1 ) + a 2 + d)d. (45) 
The numbers ao, oi, a 2 are given by 

a = J dk ( 70 Q ) 2 P+, (46a) 
ai = J dk sytf % = J dkrHiSp-Z-lEt, (46b) 

a 2 = J dk^V% = J dk( 73 Q ) 2 U 11 . (46c) 

The dagger at the column vectors denote the transposed 
vectors, i.e., row vectors. Solving Eq. 1|38[) we obtain from 



Eq. gSJ 

Q = z(a + 2a 1 )+a 2 + $(l-X)- 1 d. (47) 

The frequency dependence of Q is due to the functions F 
and K which also enter V and thus d, X and the vertex. 

A. Conductivity at zero and at high frequencies 

In the limit to — > 0, the conductivity may be deter- 
mined from the relation 

a{w = 0)=T{dQ/dz)\ z = 1 . (48) 
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Writing a 2 = za,2i + &22 and d = zD± + D2 we obtain 
with A = a + 2a! + a 2 i, 



<j(uj = 0) = t(A + D\c + ^c%=i, 



(49) 



where <? = (dc/dz) Differentiating Eq.JSSJ) with respect 
to z, 



(1 -X)c = d, 
and writing X = zX\ + X2 we obtain 
(1 - Xyg =bx + X{c. 



(50) 



(51) 



Multiplying Eg . (|50fl from the left with and using 
Eo. (|5TJ) we find for the last term in Ea. (|^|l . 

d)c = D\c + c^x\c. (52) 

Using the abbreviation cq = c\ z =i we have 

(t(u = 0)=t{A + 2D\c q +'$ X 1 c q ). (53) 

The last expression is not very transparent, but one 
can check that it exactly corresponds to the naive sym- 
metrization of 141|l using the static mean-field vertex 

r Q (k,o), 



uj 1 Im / dk / dx\dx2 



n F (xi) - n F (x 2 ) 
' uj — x\ + X2 + iO 
xTr[f Q (k, 0)GZ(xi)T a (k, 0)G£(a; a )], (54) 



in the limit uj —> 0. This result for a simpler case without 
DDW order is usually attributed to Langeri^ 

Consider now the "tail" of the Drude peak which cor- 
responds to lot 3> 1 or \z\ <C 1 and K = K(0). From 
Eq. (|45|l we obtain in linear order in z, 

Q = z(A + D\(l-X 2 )- 1 D 2 

+ d|(i-1 2 )- 1 1i(i-x 2 )- 1 5 2 (55) 

Using c\ z= o = (1 — X 2 )~ 1 D 2 we find, 

ff(wr»l) ~ -^-(A + 2D\c + ^X 1 c)\ z=0 . (56) 

UJ A T 



Eas. l4Ua|l and (|40bl) show that the two last terms in 
Ea. (|56H arc roughly by a factor JA^/t 3 smaller than the 
first three terms. Since in practice JA^/t 3 < 1 we may 
drop the last two terms in Ea. l(5Hjl and obtain approxi- 
mately, 



o{ujt > 1) ~ -5-. 

UJ^T 



(57) 



Eq. H57|) shows that the current vertex changes from its 
mean- field value Eq. I|18|) near the maximum of the Drude 
peak to its bare value in the tail of the Drude peak. The 
renormalization of the current vertex is important at low 



frequencies where it contributes to satisfy the continuity 
equation of charge. At higher frequencies this renormal- 
ization is stripped off and the bare value of the vertex 
determines the behavior in the Drude tail. This observa- 
tion is not surprizing, because the overall prefactor r _1 
shows that impurity scattering contributes in the low- 
est order. The inclusion of the scattering process into 
the corresponding current loop effectively breaks the di- 
agram into two unconnected parts with their vertices. It 
is easy to understand that these vertices are described 
by i) small but finite incoming frequency uj, ii) the full 
inclusion of the J term, and iii) the Green functions with- 
out damping r _1 . The result of the vertex summation in 
this limit is given by solving (|38[1 at z — which leads to 
T Q ~ 7" in leading order in Aq/<. 



B. Weight of the Drude peak 

The weight of the Drude peak is given by integration of 
the above a over uj for K(uj) = K(0). We define the spec- 
tral weight of the Drude peak as Wd = tt -1 / duiImQ/uj. 
The function Q(z) is an analytical function of the com- 
plex variable z = 1/(1 — iujt) in the circle \z — 1/2| < 1/2, 
which corresponds to analyticity in the upper semiplanc 
of uj. Performing a Taylor expansion around z = (a 
more accurate expansion around z = 1/2 leads to the 
same result), we write Q(z) = J2n>o QnZ n and obtain, 



Wd = I f ^L Im y 



= ^Q n = Q{z=l)-Q{z = Q). (58) 

n>l 

From Ea . (|45|l we have 

W D = A + R, (59) 
R = 0\+Dl){l-X 1 -X2)- 1 {D 1 +D 2 ) 

-Dt(l-X 2 )- 1 D 2 . (60) 

Inserting the identity 
(l-Xi-A^)- 1 = (I-X2)- 1 

+{\-X2)- 1 x 1 {\-x 1 -x 2 )-\ 

into the first term of R and using 

{i-x 1 -x 2 )- 1 d = d Q , (61) 

(l-X 2 )- 1 d = d -(l-X 2 )- 1 X 1 d , (62) 

we obtain 

R = D\c Q + Dl{l-X2)- 1 {D 1 +X l c ). (63) 

Writing D\ = d J — D\ in the second term on the right- 
hand side of Ea. l|t)3|) . using the transposed equation of 
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Ea. (|62|l . and inserting the result into Ea. l|59|l one finds 

W D = A + 2D\c q + cJXlCo 

- (c^Xt + DlXl-Xi^iDt+^co). (64) 

Without the last term in Ea. (|64|) Wd would be equal to 
the spectral weight of a simple Lorentzian of the form 
a D (u)) = <t(0)/(1 + w 2 t 2 ) with o-(Q) given by Eq.Q. 
From Ea. l|39(l follows that 1 — is a positive definite 
matrix for Ao/t <C 1. Stability arguments suggest that 
this should be generally true as long as the DDW state is 
thermodynamically stable. It follows that the last term in 
Eg. (|64fl is negative and thus reduces the spectral weight 
compared to that of a Lorentzian form for <j{lo). 

C. Optical transitions 

At finite frequencies of the order of Ao the Drudc con- 
tribution to <j(lo) is small and can be neglected whereas 
the frequency dependence of K(uj) becomes important. 
Putting z = and showing the frequency dependencies 
explicitly Eas. l|38[) and l|45|l take the form, 

c = D 2 {w)+X 2 {u)c (65) 
Q = a 22 (uj) + D\{oj)c. (66) 

If X 2 ,d 2 and 022 are purely real, then cr(uj) = ImQ/u) = 
0. This situation happens at T = 0, when to lies within 
the optical gap of K(lo). At frequencies above the thresh- 
old value, which is twice the value of the gap at the hot 
spots, i.e., to > 2A; ls , the optical conductivity is finite. 
From Eas. l|f)5|) and follows that 

Q = a 22 {cu) + D$(u)(l - X 2 (u>))- 1 D 2 (ui). (67) 

The second term on the right-hand side of Eq. (|67|l is due 
to VC and has been neglected in a previous workiAi 

D. Restricted optical sum 

Finally we discuss the restricted optical sum, given 
by the expression W = ir^ 1 f_ duiImQ/uj. From the 
Kramers-Kronig relations we find 

W = Q{lo = 0) - Q(lu = 00). (68) 

The substraction of the term Q(io — 00) reflects the fact 
that the Kramers-Kronig relation recovers the real part 
of the function up to a constant, which is Q(lu = 00). 
The latter quantity is given by taking the limit lo — > 00 
in Ea. <|67(l . From the explicit expressions Eas. (|39[l - I|40b|) 
one finds Q{to = 00) = 0. Furthermore we have, 

Q{lu = 0) = Q(z = 1) = (A + a 22 + dc)\ z=1 . (69) 

Performing a similar calculation as in the case of the 
Drude weight we find 

W = W D + (a 22 + - X 2 )- 1 D 2 )U = i. (70) 



Comparing with Ea. (|67|l one recognizes that the sum of 
the second and third terms on the right-hand side of 
Eci. (|7U|) represents just the spectral weight due to optical 
transitions as it should be. 



V. NUMERICAL RESULTS AND 
CONCLUSIONS 

Introducing the unperturbed Hamiltonian by 

H * -{ ej ' ( } 

the self-energy Sk can be written as 

After dropping the derivative V£, Ea. (|17|l represents a 
self-consistent matrix equation for Sk- The resulting gap 
equation reads, using 2 J as the energy unit, 

1 = -2 / dk 7 g(k)" F(£lk) ' nf(£2k) . (73) 

The renormalization of the diagonal part of Sk becomes 
<5£k = and \ writing <5^ k = try s (k), 7^ = ±(cos(k x ) + 
cos(ky)), 

b = 2 fdk ^(kKk ^^'^^ - (74) 

We have solved Eas. l|73|l and (|74|l by iteration. For 
T = 0,t = 0.25, t' = 0.075 we find that A is non-zero 
in the interval 0.35 < n a < 0.50, increases linearly in 
this interval from 0.15 to 0.19, and shows a step-like be- 
havior at the ends of this interval. n a is the occupation 
of one site for one spin direction, b is practically con- 
stant in the above interval and equal to —0.39. Fig|3 
shows A and 6 as a function of temperature for a fixed 
occupancy n a = 0.40. b is again practically a constant 
whereas A exhibits a BCS-like temperature dependence 
with T c - 0.07. 

Using the above values for A and b FigOshows the op- 
tical conductivity cr(w) calculated with (solid line) and 
without (dashed line) VC for the parameters indicated 
in the figure. "Without VC " means the use of the bare 
vertex, i.e., to put c 2 and C3 in Eas. (|31|) and l|32p to 
zero. In the presence of a finite A the Fermi lines are 
folded back into the small magnetic Brillouin zone, a gap 
opens near the new phase boundary in k space, and the 
resulting Fermi lines consist of arcs around the nodal di- 
rection and small pieces near the points X = (n, 0) and 
Y = (0,7r). The low-frequency Drude-like behavior of 
(t(lu) is due to electrons near the remnant Fermi lines, 
especially near the arcs around the diagonal directions^ 
The hump at frequencies above 0.2 is caused by verti- 
cal transitions between the two back-folded bands in the 
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0.00 0.02 0.04 _ 0.06 0.08 0.00 0.02 0.04 0.06 

temperature T frequency 



FIG. 2: Spectral gap Ao and the parameter & as a function 
of temperature T, using 2J as the energy unit. 



FIG. 4: Low-frequency behavior of <j(uj) with (solid curve) 
and without (dotted curve) vertex corrections. The dashed 
curve is a Lorentzian with the same value er(0) as the solid 
curve. 



t=0.25 

t'=0.075 

T=0.001 

1/x=0.01 

n o =0.40 

with vertex corr. 

without vertex corr. " 




0.0 0.1 0.2 0.3 0.4 0.5 

frequency 



FIG. 3: Optical conductivity <j(ui) with (solid line) and with- 
out (dashed line) vertex corrections. 



magnetic Brillouin zone. These optical transitions set in 
near 2A and form a broad and extended continuum due 
to transitions also far away from the Fermi lines. 

In the optical region VC increase <r(u>) slightly at fre- 
quencies not too far away from 2 A. This increase agrees 
with Eq. l|t)7|l because the second term on the right-hand 
side is due to VC and it is positive because the elements of 
X2 are small compared to unity. The influence of VC at 
small frequencies is illustrated in greater detail in FigQ] 
First, VC increase the static value cr(0) substantially. 
One reason for this is that the continuity equation for 
the charge density yields in the presence of a momentum- 
dependent order parameter an additional contribution to 
the current, given by the second term in Eq. (J2J) ■ How- 
ever, cr(0) is not only below but also above T c enhanced 
by VC. This increase is due to the renormalization of 



the one-particle energies. The resulting increase of <r(0) 
has in our approach to be ascribed to VC because only a 
non-trivial vertex allows to fulfill the continuity equation, 
see Ea. ljTH)) . Secondly, the dashed curve in Fig0] repre- 
senting a Lorentzian with the same cr(0) as for the solid 
line, illustrates the narrowing of the Drude peak with 
increasing frequency. This effect is induced by the fre- 
quency dependence of the vertex. The figure shows that 
the solid line approaches at larger frequencies the curve 
without VC (dotted line) which is a Lorentzian. This 
crossover takes place at the impurity scattering rate 1/r. 
FigEl a l so shows that the effect of VC cannot be simu- 
lated in a theory without VC by a suitable redefinition 
of the current. The use of the bare current (first term 
in Eq.||2J)) would result in a violation of the continuity 
equation and miss the enhancement of cr(w) at low fre- 
quencies. On the other hand, the use of a current derived 
from the continuity equation as in Eq. would yield the 
dashed curve in Fig0]and would thus overestimate <j{oj) 
at large frequencies. 

FigEl shows the temperature dependence of the re- 
stricted optical sum W and its low-frequency part, the 
Drude weight Wjj. VC affect W in a two- fold way: They 
enhance the weak temperature dependence of W, cal- 
culated without VC (dot-dashed line), and also increase 
substantially the absolute value of W. The optical spec- 
tral weight W — Wd is expected to set in below T c in a 
square-root manner and to be rather independent of VC. 
Thus the enhancement of W and its stronger temperature 
dependence by VC should be related to the low-frequency 
part of <t(w), i.e., to Wd, which is in agreement with the 
dashed and dotted curves in Fig[SJ 

Decreasing t' means that the length of the remnant 
Fermi lines becomes smaller until they shrink to points 
along the diagonals for t' = as in a d-wave superconduc- 
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FIG. 5: Restricted optical sums with (solid line) and without 
(dot-dashed line) vertex corrections, and Drude weights with 
(dashed line) and without (dotted line) vertex corrections. 




0.2 0.3 

frequency 



FIG. 6: Optical conductivity a(cu) for T = 0, A = 0.197 (solid 
line), T = 0.06, A = 0.158 (dashed line), and T = 0.08, A = 
0.041 (dot-dashed line). 



only a small part of W is left for the Drude peak. With 
increasing temperature the length of the Fermi lines in- 
creases and the spectral weight is shifted from the optical 
region to the Drude peak, as shown by the dashed and 
dot-dashed lines in this figure. 

In conclusion, we have calculated the optical conduc- 
tivity <j{u>) of a two-dimensional d-CDW conductor tak- 
ing into account vertex corrections. The employed for- 
mulation can easily be generalized to an arbitrary sym- 
metry of the CDW state and to a general band structure. 
We have addressed and solved the problem of a proper 
definition of the current to be used in calculating <j(u>). 
Defining the current via the continuity equation both the 
renormalization of the one-particle energies due to the in- 
teraction causing the CDW as well as an additional term 
due to the momentum dependence of the CDW order 
parameter enter. Using this renormalized current and 
bare vertices yields a o~(w) which is shown to be correct 
only in the limit ix> — > 0. On the other hand, using bare 
vertices and the bare electron dispersion violates in gen- 
eral the continuity equation for the charge density and 
gives a <j(ui) which is correct only at large frequencies. 
The proposed solution is to use bare current vertices but 
to include vertex and the corresponding self-energy cor- 
rections. In this way the continuity equation is explicitly 
fulfilled and the frequency dependence of the vertex inter- 
polates smoothly from the strongly renormalized regime 
at low to the unrenormalizcd regime at high frequencies. 
This behavior is also visible in the numerically calcu- 
lated conductivity curves. These calculations show the 
spectral features of the Drude-like peak at low and the 
transitions across the gap at high frequencies and their 
dependence on temperature and the next-nearest neigh- 
bor hopping constant t'. Remarkable is also that the 
temperature dependence of the restricted optical sum W 
(i.e., the spectral weight integrated over all frequencies) 
between T = and T = T c is enhanced by vertex cor- 
rections. This may be interpreted as a violation of the 
sum rule by about 10 per cents in that interval. This as 
well as the other above predictions may help to identify 
experimentally unconventional CDW's. 



tor. Thus the ratio Wd/W should decrease with t'. The 
solid line in FigHj] calculated with t' = 0.05, T = 0, n a — 
0.48, 1/r = 0.03, illustrates that this indeed is the case. 
Most of the spectral weight resides now in a very broad 
and extended high-frequency structure above 2A ~ 0.35, 
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